ON COLLISIONS OF BROWNIAN PARTICLES 



TOMOYUKI ICHIBA AND IOANNIS KARATZAS 

Abstract. We examine the behavior of n Brownian particles diffusing 
on the real line with bounded, measurable drift and bounded, piecewise 
continuous diffusion coefficients that depend on the current configura- 
tion of particles. Sufficient conditions are established for the absence, 
and for the presence, of triple collisions among the particles. As an 
application to the Atlas model for equity markets, we study a special 
construction of such systems of diffusing particles using Brownian mo- 
tions with reflection on polyhedral domains. 



1. Introduction 

It is well known that, with probability one, the n-dimensional Brown- 
ian motion started away from the origin will hit the origin infinitely often 
for n = 1, while it will never hit the origin for n > 2. This is also true 
for the n— dimensional Brownian motion with constant drift and diffusion 
coefficients, by Girsanov's theorem and re-orientation of coordinates. The 
next step of generalization is the case of bounded drift and diffusion coef- 
ficients. The existence of weak solutions for the stochastic equations that 
describe such processes was discussed by Stroock & Varadhan [24], Krylov 
|17j through the study of appropriate martingale problems. 

Now let us suppose that W 1 is partitioned as a finite union of disjoint 
polyhedra. Bass & Pardoux [3] established the existence and uniqueness of 
a weak solution to the stochastic integral equation 

(1.1) X(t)=x +f v(X(s))ds+ [ a(X(s))dW(s), < t < oo 
Jo Jo 
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with initial condition xq G M n , where the measurable functions \x : W 1 — ► IR n 
and cr : M n — > M nxn are bounded and, moreover, a is everywhere non- 
singular and piecewise constant (that is, constant on each polyhedron). The 
continuous process {W(i), < t < 00} is an n— dimensional Brownian mo- 
tion on some filtered probability space (Q, {Ft}-, IP). Here uniqueness is 
understood in the sense of the probability distribution. 

Bass and Pardoux also discovered an interesting phenomenon, namely, 
that the weak solution to (11. ip may satisfy 

(1.2) F Xo (X(t) = 0, i.o.)=l; xq g M n , 

for a diffusion matrix cr(-) with special structure and without drift //(•) = . 
Here ¥ Xo is the solution to the martingale problem corresponding to (jl.ip . 
In the Bass &: Pardoux [3] example the whole space W 1 is partitioned into 
a finite number of polyhedral domains with common vertex at the origin, 
carefully chosen small apertures, and cr(-) constant in each domain. We 
review this example in Remark 12.41 

In the present paper we find conditions sufficient for ruling (|1.2p out. 
More specifically, we are interested in the case of a bounded, measurable 
drift vector //(•) and of a bounded, piecewise continuous diffusion matrix 

m 

(1.3) a{x) = ^o u (x)l nv (x) = a p[x) {x)- x G M. n , 

v=l 

under the assumption of well-posedness (existence and uniqueness of solu- 
tion) when n > 3. Here is the indicator function; the sets {TZ^}™^ 
form a partition of R n for some m G N , namely, 1Z U n 1Z K = for v ^ k 
and U^ =1 TZ U = M. n ; and the mapping p : M. n — ► {1, . . . , m} satisfies x G 7Zp(x) 
for every x G W 1 . Throughout this paper we shall assume that 1Z U is 
an n-dimensional polyhedron for each v = 1, . . . , m , and that the (re x re) 
matrix-valued functions {^(•) "i/(')} y= i are positive-definite everywhere. 

We shall also assume throughout that there exists a unique weak solution 
for the equation (II. ip . Existence is guaranteed by the measurability and 
boundedness of the functions /x(-) and <r(- )(/(•) , as well as the uniform 
strong non-degeneracy of a{-)a'{-) (eg. Krylov [16], Remark 2.1), where the 
superscript / represents the transposition. Uniqueness holds when re = 1 or 
n = 2 ; for n > 3 , the argument of Chapter 7 of Stroock & Varadhan [21] 
implies uniqueness if the function er(-) in (II. 3D is continuous on M n (The- 
orem 7.2.1 of |24|) or close to constant (Theorem 7.1.6 of [24]), namely, if 
there exists a constant (n x n) matrix a and a sufficiently small 5 > 
depending on the dimension n and the bounds of eigenvalues of cr(-) such 
that \\a(x)a'(x) — a\\ < 5 . Bass Sz Pardoux [3] showed uniqueness 

for piecewise-constant coefficients, i.e., ov(-) = u v , v = 1, . . . , m. For fur- 
ther discussion on uniqueness and non-uniqueness, we refer to the paper by 
Krylov [16] and the references in it. The structural assumption (jl.3p may 
be weakened to more general bounded cases, under modified conditions. 
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Our main concern is to obtain sufficient conditions on /i(-), and on <r(-) 
of the form f)l .3f) . so that with n > 3 we have 

P xo (Xi(t) = Xj{t) = X k (t) , for some t > 0) = or 

lLI ' P x . (Xi(t) = Xj(t) = X k (t) , for some t > 0) = 1 ; xq G R n 

for some l<i<j<k<n. Put differently, we study conditions on 
their drift and diffusion coefficients, under which three Brownian particles 
moving on the real line can collide at the same time, and conditions under 
which such "triple collisions" can never occur. Propositions 1 and 2 provide 
answers to these questions, in Section El 

In Section [3] we study a class of the weak solutions to the stochastic differ- 
ential equation (jl.lj) . clarifying the relationship between the rank of process 
coordinates and the reflected Brownian motion on (n — 1)— dimensional poly- 
hedral domain. Proposition 3 shows that the process has no triple collisions 
under some parametric conditions. 

The results have consequences in the computations of local times for the 
differences {Xi(t) — Xj{t),Xj{t) — Xk(t)}. We discuss such local times with 
application to the analysis of a so-called "Atlas model" for equity markets 
in Section 01 Proofs of selected results are presented in Section 

Recent work related to this problem was done by Cepa & Lepingle [5]; 
they consider a system of mutually repelling Brownian particles and show the 
absence of triple collisions. The electrostatic repulsion they consider comes 
from unbounded drift coefficients; in our setting, all drifts are bounded. 

2. A First Approach 

2.1. The Setting. Consider the stochastic integral equation (jl.ip with 
coefficients u(-) and cr(-) as in (|1.3j) . and assume that the matrix- valued 
functions c u (-), u = 1, . . . ,m are uniformly positive-definite. Then the in- 
verse cr _1 (-) of the diffusion coefficient cr(-) exists in the sense cr _1 (-) = 
Y^=\ a u 1 (')^-Tl^(') ■ As usual, a weak solution of this equation consists of a 
probability space (SI, J 7 , P) ; a filtration {Tt , < t < oo} of sub-cr-fields of T 
which satisfies the "usual conditions" of right-continuity and augmentation 
by the IP— negligible sets in T\ and two adapted, n-dimensional processes on 
this space X(-), W(-) on this space, such that W(-) is Brownian motion 
and (jl.ip is satisfied P— almost surely. The concept of uniqueness associated 
with this notion of solvability, is uniqueness in distribution for X(-) . 

2.2. Removal of Drift. We start by observing that the bounded drift has 
no effect on the probability of absence of triple collisions. Indeed, if we 
define an n— dimensional process := a^ 1 (X (t)) fj>(X (t)) , < t < oo , 
then the nature of the functions fi(-) and <r(-) in (jl.3p guarantees that the 
mapping t \— » ^(t) is right-continuous or left-continuous on each bound- 
ary &R,p(x(t)) at each time t, deterministically, according to the position 
Tlp(x{t-)) °f X(t— ). Thus, although the sample path of n-dimensional pro- 
cess £(•) is not entirely right-continuous or left-continuous, it is progressively 
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measurable. Moreover, £(•) is bounded, so the exponential process 



(2.1) r/(i) = exp 



< t < oo 



is a continuous martingale, where ||x|| 2 := X^i^f; x £ ^ n stands for n- 
dimensional Euclidean norm and the bracket (a;, y) := ^YTj=\ x jVj ^ s the 
inner product of two vectors x,y S R n . By Girsanov's theorem 

(2.2) W{t):=W{t)+ I (j- l {X(u))Li{X{u))du, T t \ < t < oo 

is an n-dimensional standard Brownian motion under the new probability 
measure Q, locally equivalent to P , that satisfies 

(2.3) Q X0 (C) = E p *o ( V (T)1 C ); Cef T , < T < oo . 

Let us define an increasing family of events Ct '■= {Xi(t) = Xj(t) = 
Xk(t) , for some t E [0, T] } , T > . If we know a priori that 

(2.4) Q X0 {Xi(t) = Xj(t) = X k {t) , for some t > 0) = , 

then we obtain = Q X() (Q) = ^x Q {Ce) for £ > 1 , and so 

Pxo = = X k (t) , for some t > 0) = F^U^d) 

( 2 - 5 ) =limF IO (Q)=0. 

Thus, in order to evaluate the probability of absence of triple collisions in 
(jl.4p . it is enough to consider the case of //(•) = in (jl.ip . namely 

(2.6) X(t)=x +[ a(X(s)) dW(s) , < t < oo , 

■/ o 

under the new probability measure (Q^q • The infinitesimal generator A of 
this process, defined on the space C 2 (R n ;R) of twice continuously differen- 
tiable functions ip : R n —> R , is given as 

1 n (9 2 

(2.7) Aip(x) :=-^ OifcCx) -Q-^T k M*)} \ V G ^ 2 (R n ; R), 

i,/c=l 

where fjj(-) is the (i, j)-th element of the matrix-valued function cr(-) and 

n 

(2.8) a ifc (x) := ^a^x) a kj {x) , A(x) := {a^)}i<ij< n 5 x G M ™- 

The uniform positive-definiteness of the matrices {(t^(t^}(-), f = 1, • • ■ , m in 
(jl.3p implies that the operator A is uniformly elliptic. As is well known from 
[24] . existence (respectively, uniqueness) of a weak solution to the stochastic 
integral equation (|2.6p . is equivalent to the solvability (respectively, well- 
posedness) of the martingale problem associated with the operator A . 
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2.3. Comparison with Bessel processes. Without loss of generality we 
start from the case i = l,j = 2,k = 3'm (|1.4p . Let us define (n x 1) 
vectors di,d 2 ,d 3 to extract the information of the diffusion matrix cr(-) on 
(X 1 ,X 2 ,X 3 ), namely 

d l := (1, -1, 0, • • • , 0)', d 2 := (0, 1, -1, 0, • • • , 0)', d 3 : = (-1, 0, 1, 0, • • • , 0)', 

where the superscript / stands for transposition. Define the (n x 3)-matrix 
D = (d±, d 2 , d 3 ) for notational simplicity. The cases we consider in (jl.4p for 
i = 1, j = 2, k = 3 are equivalent to 



F Xo [s 2 (X(t)) = , for some t > 0J = and 
P X0 (s 2 (X(t)) = , for some t > o) = 1 ; xq £ 



where the continuous function 

s 2 (x) := (xi - x 2 ) 2 + (x 2 - x 3 ) 2 + (x 3 - xi) 2 

= d[xx'di + d' 2 xx'd 2 + d' 3 xx'd 3 = x'DD'x ; x € M™ 

measures the sum of squared distances for the three particles of interest. 
Thus, it suffices to study the behavior of the continuous, non-negative pro- 
cess {s 2 (X(t)); < t < oo} around its zero set 

(2.10) Z := {x G R n : s(x) = 0}. 

Let us define the following positive, piecewise continuous functions Q(-) , 
R(-) computed from the variance-covariance matrix A(-) = a(-)a'(-) : 

~ _ tr&ce(D'A(x)D) ■ x'DD'x _ trace(D'A(x)D) 
(2n) {X}: ~ x'DD'A(x)DD'x " Q{xj ' 6 
^ x'DD'A{x)DD'x mn . „ 

Under the new probability measure Q xo of (|2.3|) the process s(X(-)) is a 
semimartingale with decomposition ds(X(t)) = h(X(t))dt + d@(t) , where 

1 ( 3 II 3 

M^) : = 2 3 / \ s2 (^) Yl d ' i a{x)a(x)' di -\\Y^ a ( x Y d i4" 

s ^ V i=i i=i 

/ 2 12 n _ x'DD'x ■ trace (Z?'i4(a;).D) - x'DD'A(x)DD'x 

~ 2{x 7 DD i xjy~ 2 



©(*) 
<©>(*) 



(fi(x) -l)Q(x) 
2s(x) 

S (X(r)) 
* x'DD'A(x)DD'x 




o 



x'DD'x 




and 



dr= / Q(X(r))dr; < t < oo . 

s=X(r) 
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respectively. Here, as we shall see (I2.27P in Remark 12. 1\ 

(2.13) Q(-)>c :=3 min A,(x) > in R n \ Z 

l<i<n,xm n \Z 

for the eigenvalues {Aj(-),1 < i < n} of ^4(-) , and so (0)(-) is strictly- 
increasing, when X{-) G W 1 \ Z . Now we define the increasing family of 
stopping times A u := inf{t > : (0) (t) > u} , < u < oo and note that 
we have 

s(u) := s(X(Au)) = s(x ) + / h(X(t))dt + B(u); < u < oo , 

Jo 

where B(u) := Q(A U ) , < u < oo is a standard Brownian motion, by the 
Dambis-Dubins-Schwartz theorem on time-change for martingales. Thus, 
with d(u) := R(X(A U )) we can write 

(2.14) ds(u)= ^l^ 1 du + dB{u); < u < oo , 
because 

[R(X(Au)) - 1] Q(X(Au)) 1 5(«)-l 



/ l (X(A u ))A; - 2s (X(A u )) Q(X(A n )) - 2*(u) • 

The dynamics of the process s(-) are therefore comparable to those of the 
5— dimensional Bessel process, namely 

dt(u) = - \ du + dB{u) ; < u < oo . 
2 t(u) 

By a comparison argument similar to Ikeda &; Watanabe [14] and Exercise 
5.2.19 in [15] , we prove in subsection 15.11 the following result. 

Lemma 2.1. Suppose xq G W 1 \ Z . If d := essinf info<t <00 d(t) > 2 , 

(2.15) Qj;,, ( s(t) > , for some t > ) = . 
If, on the other hand, d := essup sup 0<i<oo d(t) < 2, then 

(2.16) Q xo ( s(t) = , for infinitely many t > ) = 1 ; 
and we /iaue t/ie following estimate 

(2.17) Q Xo (s(t) = 0, for some t G [0,T]) > 1 - k(T;s(x ),0) 

where n(- ; y, 5) is i/ie tail distribution of the first hitting-time at the origin 
for Bessel process in dimension 5 G (0, 2) , starting at y > : 

I 00 1 / v 2 \ 5 v 2 

(2.18) K(T',y,S):=J Jf{5)\2t) dt; ° " T < °° ' y > °' 



This function decreases as T~ s with T f oo. Combining Lemma 12.1 
with the reasoning in subsection 12.21 and the definition d(-) = R(X(A.)) , we 
obtain the following result on the absence of triple collisions: 
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Proposition 1. Suppose that the matrices cr v ('), u ,m in (ji.ffl are 

uniformly bounded and positive-definite, and satisfy the following condition 

(2.19) inf R(x) > 2 

xm n \z 

for R(-) in \2.1 1\ . Then for the weak solution X(-) to (\2.6\i we have 

Q XQ = X 2 (t) = X 3 (t) , for some t > 0) = , V xo £ R" \ 2 . 

Reasoning as in pO )- p£S| /or t/ie u>ea£; solution X(-) to HI. 3 . we get 

(2.20) P^ (Xi(i) = X 2 (i) = X 3 (f) , for some t > 0) = , V x e E n \ Z . 

A class of examples satisfying (|2.19p is given in Remarks I2.2H2.3I and 
Section 15.31 below. On the other hand, regarding the presence of triple 
collisions, we have the following result; its proof is in subsection 15.21 

Proposition 2. Suppose that the matrices o~ u (-), v = 1, ■ • ■ ,m in d~OJ) are 

uniformly bounded and positive-definite, and 

(2.21) 5 : = sup R(x) < 2 . 

x£R"\Z 

Then the weak solution X(-) to (|j?.6]) starting at any xq € M. n satisfies 
Qx (*i(t) = X 2 (t) = X 3 (t), for some t > 0) = 1 , 

and we have an estimate similar to f\2.17\) : 
(2.22) 

Q X0 (X 1 (t) =X 2 (t) = X 3 (t), for some t€ [0, T]) > 1 - k(c T ; s(x ), <5 ) • 

/fere i/ie distance function s(-) and i/ie tail probability K(-;y,5o) are given 
by (|j?..9p and i<Sj) . now lozi/i dimension 5o £ (0,2) as in \2.21\ . and the 
positive constant cq is given by <\2.13^ . 

Moreover, if <5* := supj. 6B nw R(x) < 2 holds for the modification 

[tmce(D'A(x)D) + 2x'DD'fi(x)} ■ x'DD'x 



R(x 
(2.23) 



x'DD'A(x)DD'x 



of the function R(-) in (\2.11\) , then 

(2.24) P X0 (A"i(t) = X 2 (t) = X 3 (t), for some t > 0) = 1 

and we have an estimate similar to \2.1T§ , \2.22\ : 
(2.25) 

P X0 ( X x (t) = X 2 (t) = X 3 (t), for some t € [0, T]) > 1 - k(c T ; s(x ), <5*) • 

Remark 2.1. Since A(-) is positive-definite and rank(D) = 2, the matrix 
D'A(-)D is non-negative-definite and the number of its non-zero eigenvalues 
is equal to rank(D'A(-)D) = 2 . This implies 

ELaP(x) 

maxi<j< 3 Af (x) 



R(x) > ^ i=1 ' V ^ - > 1; iel"\2, 
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where {Af (•),« = 1,2,3} are the eigenvalues of the (3 x 3) matrix D'A(-)D. 
On the other hand, an upper bound for R(-) is given by 



(2.26) 



R(x) < 



trace (D'A(x)D) 
3mini<j< n Aj(x) 



x G 



where {Aj(-), 1 < i < n} are the eigenvalues of A(-). In fact, we can verify 
DD'DD' = 3DD', {x G W 1 : DD'x = 0} = Z, and so if DD'x / £ W 1 , we 
obtain the upper bound (|2.26j) for R(-) from 

. w w x'DD'A{x)DD'x Q(x) trace (D'i(x)D) 

(2.27) mm Ai(x) < — — — — — = — - — = . □ 

v ; i<i<n n ; - x'DD'DD'x 3 3_R( X ) 

Remark 2.2. For the standard, n— dimensional Brownian motion, i.e., cr(-) = 
J n , n > 3 , the quantity R{-) of (|2.1ip is computed easily: R(-) = 2 . More 
generally, suppose that the variance covariance rate A(-) is 



A(x) := ( a ^n + P*DD' + II' diag( 7l ,)) • l Uu (x) ; xGK", 
i/=i 

for some scalar constants a u , (3 V and {n x 1) constant vectors j u , f = 
1, . . . , m. Here diag(x) is the (n x n) diagonal matrix whose diagonal entries 
are the elements of x G M n , and I is the (n x 1) vector with all entries equal 
to one. Then R(-) = 2 in M n \ Z , because I'D = (0, 0, 0) G M lx3 and 



= -DD'DD' 
3 



V 



-1 -1 

2 -1 

-1 2 








Hence, if the coefficients a u ,f3 u and j v , v = 1, . . . ,m are chosen above so 
that A(-) is positive-definite, we have (|2.2U|) . □ 



Remark 2.3. The condition (|2,19[) in Proposition 1 holds under several cir- 
cumstances. For example, take n = 3 and fix the elements an(-) = a 22 (-) = 
a 33 (-) = 1 of the symmetric matrix ^4(-) = crcr'(-) in (|2.8|) and choose the 
other parameters by 



a 12 (x) = a 2 i(x) 
(2.28) a 23 (x) = a 32 (x) 
a 31 (x) = a 13 (x) 



: a 2+ ln 2+ (x) + a 2 -l^ 2 _ (x) , 

: "3+ 1^3+ (x) + a 3 _ 1^ 3 _ (x) ; x G 



where 


1Zi±, i 


= 1 




= { x G 


R 3 


Ki- 


= {ie 


R 3 




= {ie 


R 3 


^3- 


= { x G 


M 3 


h(x) 


= [X3 ~ 


Xl 




= [X2 ~ 


X3 


h(x) 


= [x\ - 


X2 
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!, 3 are subsets of R 3 denned by 

h(x)>0}, U 2+ :={x£R 3 : fi(x) = , f 2 (x) > } 
fl(x)<0}, ft 2 _ :={xGM 3 : fr(x) = , f 2 (x) < } 
fi(x) = f 2 (x) = 0,f 3 (x) >0}, 
fl(x) = f 2 (x) = 0,f 3 (x) <0}, 
(-2 + \/3)(x 2 - x 3 )] • [x 3 - xi 



Xl 


-(-2- 


V3)(x 2 


- »3)] , 


X3 


-(-2- 


73) (xi 




X 2 


-(-2- 


V3)(x 3 




tig 


< a i+ 


< 1/2, 


-1/2 < 



for x G R with 1 

«i_ < for % = 1,2,3. Note that the zero set Z defined in (I2.10P is 
{ x G R 3 : fi(x) = fz(x) = f 3 (x) = } . Thus, we split the region R 3 \ Z into 
six disjoint polyhedral regions 1Zi± ,i = 1, 2, 3 . See Figure EJ and subsection 
5.31 for the details of this example. □ 

Remark 2.4. In the example of Bass & Pardoux [3], mentioned briefly in 
the Introduction, the diffusion matrix a(-) = Y1T=1 (') i n ([HH]) has 

a special characteristic in the allocation of its eigenvalues: All eigenvalues 
but the largest are small, namely, they are of the form (l,e, • • • ,e), where 
< e < 1/2 satisfies, for some < 5 < 1/2 : 

x / cj(x)(t'(x)x 



Ixll 2 



< 5 for x G M n \ {0} , and ^ ^ + - < 1 

l — o 



This is the case when the diffusion matrix er(-) can be written as a piecewise 
constant function X)I/Li IJ v^-'R.v(') ; where the constant (n x n) matrices 
{ov, i/ = 1 . . .m} have the decomposition 



g v o' v := (y u , B v ) diag (l, e 2 , . . . , e 2 ) ( 
the fixed (n x 1) vector G Mj, satisfies 

hu\\ = l, 1{X ' f 2 )|2 >l-g; xg^\W, 
||x|| z 

and the (nx(n-l)) matrix £>j, consists of (n— 1) orthonormal n— dimensional 
vectors orthogonal to each other and orthogonal to y v , for v = 1, ...m. 
Then, for all x G M n , we have 

llxll 2 trace (a(x)a'(x)) , (n — l)e 2 + (5 

V V / V H Y < ^ 1 < 1. 



x'cr^cr'^x 1 — 5 

This is sufficient for the process X(-) to hit the origin in finite time. 

To exclude this situation, we introduce the effective dimension ED^-) of 
the elliptic second-order operator A defined in (|2.7p . namely 

||x|| 2 trace (a(x)a' (x)) ||x|| 2 trace (A(x)) 



(2.29) ED4(x) 



x'a(x)a'(x)x x'A(x)x 
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for x S W 1 \ {0}. This function comes from the theory of the so-called exte- 
rior Dirichlet problem for second-order elliptic partial differential equations, 
pioneered by Meyers & Serrin |18j . These authors showed that 

(2.30) inf EBa(x) > 2 

xeK»\{o} 

is a sufficient condition for the existence of solution to an exterior Dirichlet 
problem. In a manner similar to the proof of Proposition 1, it is possible to 
show that (j2T30j) is sufficient for F Xo (Xi(t) = ■ ■ ■ = X n (t) = for some t > 
0) = , since R(-) becomes ED^(-) when the matrix D is replaced by the 
identity matrix. (In this manner, the function R(-) of (|2.12p is interpreted 
as a "local" version of the effective dimension.) 

With a(-) as in (jl.3p . the effective dimension ED^(-) satisfies 

t-it-v , > . • f \\x\\ 2 trace {a v {x)a' u {x)) \ . ( £" =1 Ai„(a;) 

hU^(x) > mm I — > mm 1 



=1,— ,m\ x' a v {x)a' v {x)x J v=l,— ,m ymaxj^... 5 „ Xi u (x] 

for x £ W 1 \ {0}, where {Ai^(-), i = 1, ••• ,n} are the eigenvalues of the 
matrix-valued functions o~ u (-)o~' u (-), for v = 1, ■ ■ ■ , m . Thus, ED^(-) > 2 if 

• r ( ELl-M^) \ ^ o 

mi mm — — -— > 2 ; 

xeK"\{o} v=l,-,m ^max^i^..^ Xi u {x) ) 

this can be interpreted as mandating that the relative size of the maximum 
eigenvalue is not too large, when compared to all the other eigenvalues. □ 

Remark 2.5. Friedman [9] established theorems on the non-attainability of 
lower-dimensional manifolds by non-degenerate diffusions. Let M be a 
closed k— dimensional C 2 — manifold in M n , with k < n — 1 . At each point 
x £ M , let Afc + j(a;) form a set of linearly independent vectors in R" which 
are normal to M at x . Consider the matrix a{x) := (otij(x)) with 



n 



ij 



(x) = (A(x)N k+i (x),N k+j (x)) ; 1 < i, j < n - k , x G M 



Roughly speaking, the strong solution of (jl.ip under linear growth condition 
and Lipschitz condition on the coefficients cannot attain M , if rank (a(x)) 
> 2 holds for all x G M . The rank indicates how wide the orthogonal 
complement of M is. If the rank is large, the manifold Ai is too thin to 
be attained. The fundamental lemma there is based on the solution u(-) of 
partial differential inequality Au(-) < fJ,u(-) for some fi > , outside but near 
M with linidis^ ^m)-^oo u(x) = oo , which is different from our treatment in 
the previous sections. 

Ramasubramanian [20], [21] examined the recurrence and transience of 
projections of weak solution to (jl.ip for continuous diffusion coefficient cr(-) , 
showing that any (n — 2)— dimensional C 2 — manifold is not hit. The integral 
test developed there has an integrand similar to the effective dimension 
studied in [18] . as pointed out by M. Cranston in Mathematical Reviews. 

Propositions 1 and 2 are complementary to these previous general results, 
since the coefficients here are allowed to be piecewise continuous; however, 
they depend on the typical geometric characteristic on the manifold Z we 
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are interested in. Since the manifold of interest in this work is the zero 
set Z of the function s(-) , the projection s(X(-)) of the process and the 
corresponding effective dimensions ED^(-) and R(-) are studied. □ 

Remark 2.6. As V. Papathanakos first pointed out, the conditions (|2.19p . 
(|2.2ip in Propositions 1 and 2 are disjoint, and there is a "gray" zone of sets 
of coefficients which satisfy neither of the conditions. This is because we 
compare with Bessel processes, replacing the n-dimensional problem by a 
solvable one-dimensional problem. In order to look at a finer structure, we 
discuss a special case in the next section by reducing it to a too-dimensional 
problem. This follows a suggestion of A. Banner. □ 

3. A Second Approach 

In this section we discuss a class of weak solutions to the equation (jl.ljl 
with the structure (11.31) which exhibit "no triple collisions''' , using the n— dimensional 
ranked process and the (n — 1)— dimensional reflected Brownian motion on 
polyhedral domains. 

3.1. Ranked process. Given a vector process X(-) := {(Xi(t), . . . , X n (t)) ; 
< t < oo}, we define the vector Xt.\ := {(X(u(t), . . . ,Xt n \(t)) ;0 < t < 
oo } of ranked processes, ordered from largest to smallest, by 

(3.1) X (k) (t):= max ( min . . . , X ih {t))) ; < t < oo 

l<ii<. ..<ifc<n V / 

for k = 1, ■ ■ ■ , n . If, for every j = 1, . . . , n — 2, the two-dimensional process 

(3.2) toO, Y j+1 (-))' := {X {j) (-) - X (j+1) (-), X [J+1) (.) - X {j+2) (-))' 

obtained by looking at the "gaps" among the three adjacent ranked processes 
X(j\(-), X(j+i)(-), Xtj + 2)(~) i never reaches the corner (0,0)' of E? , almost 
surely, then the process X(-) satisfies 

(3.3) F X0 (Xi(t) = Xj(t) = X k (t) , for some k) , t > 0) = 

for xq E R n \ Z . On the other hand, if for some j = 1, . . . , n — 2 the vector 
of gaps — Xfj + i)(-), Xrj + i\(-) — Xh+2)('))' does reach the corner 

(0, 0)' of M. 2 almost surely, then we have 

F Xo (Xi(t) = Xj(t) = X k (t), for some (i, j, k), t > 0) = 1 ; x G R n . 

Thus, we are led to study the ranked process Xt.) and its adjacent differ- 
ences. In the following we use the parametric result of Varadhan & Williams 
|25j on Brownian motion in a two-dimensional wedge with oblique reflection 
at the boundary, and the result of Williams [26] on Brownian motion with 
reflection along the faces of a polyhedral domain. 

There is a long list of contributions to the study of attainability of the ori- 
gin for the Brownian motion with reflection. Recently Delarue [6] considered 
the hitting time of a corner by a reflected diffusion in the square. Rogers |22j 
[23] . Burdzy & Marshall [4] considered Brownian motion in a half-space with 
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variable angle of reflection. Here we consider oblique constant reflection on 
each face of the polyhedral region. 

3.2. Reflected Brownian motion. Let ei,...,e n _i be unit vectors in 
R n , n > 3 , and consider the non-negative orthant 

S := R^- 1 = | Y^ x k£k ■ xx > 0,... , x„_i > | 

whose (n — 2)— dimensional faces 3l, ... >3n-i are given as 

fn-l \ 

$i := < Xkek : Xk > for k = 1, . . . , n — 1, X{ = > ; 1 < i < n — 1 . 

Let us denote the (n— 3)— dimensional faces of intersection by := S^H^' 
for 1 < i < j < n — 1 and their union by ^° := ^>\<i<j<n-vB°ij ■ 

We define the (n — 1)— dimensional reflected Brownian motion Y(-) := 
{(Yi(t), . . . , ; t > 0} on the orthant M™^ 1 with zero drift, constant 

((re — 1) x (re — 1)) constant variance/covariance matrix 21 := EE', and 
reflection along the faces of the boundary along constant directions, by 

(3.4) Y(t) = y(0) + EB(t) + 5RL(t) ; < t < oo , Y(0) E M"" 1 \ £° . 

Here, {-B(t); < t < oo} is (n — 1)— dimensional standard Brownian motion 
starting at the origin of M n_1 . The ((n — l)x(re — 1)) reflection matrix 
y{ has all its diagonal elements equal to one, and spectral radius strictly 
smaller than one. Finally, the components of the (n — 1) — dimensional pro- 
cess L(t) := (Li(t), . . . ,L n _i(i)); < t < oo are adapted, non-decreasing, 
continuous and satisfy Yi(t) dLi(t) = (that is, £«(•) is flat off the set 
{t > : Yi(t) = 0}) almost surely, for each i = 1, ... re — 1. Note that, if 
Y(t) lies on = n 3j , then Yi(t) = Y s (t) = for 1 < i ^ j < n - 1. 

Harrison & Reiman jTT] introduced and constructed this process pathwise, 
through the multi-dimensional Skorohod reflection problem. 

3.2.1. Rotation and Resettling. Assume that the constant covariance matrix 
21 = ££' is positive-definite; let U be a unitary matrix whose columns 
are the orthonormal eigenvectors of 21 ; and let £ be the corresponding 
diagonal matrix of (positive) eigenvalues such that £ = f/'2t U . Define 
y(-) := £ _1 / 2 C/y(-) and note that, by this rotation and rescaling, we obtain 

Y(t) = Y (0) + B(t) + £- 1/2 UmL(t) ; < t < oo 

from (|3.4p . where 5(-) := £ _1 / 2 [/Si?(-) is another standard (n— 1)— dimensional 
Brownian motion. We may regard Y"(-) as reflected Brownian motion in a 
new state-space <S := H^^-^UMJ^r 1 . The transformed reflection matrix 
can be written as 

« = fi- 1 / 3 ^ = (^ + Q)C = (ti, . . . , tn-i) , where 
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£: = 2r 1/2 , Q : = diagffl) , 91 := £ 1/2 [/£ = (n h . . . , n n _i) , 
(3 ' 5) O := ST^l/Wr 1 - VI = (q a , . . . , q„_ a ) . 

Here 2) = diag{%) is the ((n — 1) x (n — 1)) diagonal matrix with the same 
diagonal elements as those of 21 = EE' (the variances). The constant vectors 
ti, cfi, Hj, i = 1, . . . , n — 1 are ((n — 1) x 1) column vectors. 

Since U is an orthonormal matrix that rotates the state space & = R'™ -1 , 
and £ 1//2 is a diagonal matrix which changes the scale in the positive direc- 
tion, the new state-space S is an (n — 1)— dimensional polyhedron whose 
z-th face := £~ l ' 2 U$i has dimension (n — 2), for i = 1, . . . , n — 1. 

Note that diagfii'Q) = and diag(9T'9T) = /, that is, ttj and c|j are 
orthogonal and is a unit vector, i.e., n^q, = and n^rij = 1 for i = 
1, . . . , n — 1. Also note that rij is the inward unit normal to the z-th face & 
of the new state-space S , on which the continuous, non-decreasing process 
Lj(-) actually increases, for i = 1, . . . , n — 1. The i-ih face can be written 
as {x G 6 : n^x = bj} for some bj G R, for « = 1, . . . , n — 1. 

Moreover, the i-th column tj of the new reflection matrix is decomposed 
into components that are normal and tangential to , i.e., t$ = £a(x\i + c[i) 
for i = 1, ... n — 1, where (£u is the (i, i)-element of the diagonal matrix £. 
Since the matrix ST'^'^U of the transformation is invertible, we obtain 

(3.6) ¥{■) G : = & n ^ ^ e 3y ; 1 < i < j < n - 1 . 

Thus, in order to decide whether the process Y(-) in (|3.4p attains #° , it 
is enough to decide whether the transformed process Y(-) attains the set 
5° := £ l / 2 U$° = Ul<i<i<n-1 $ij ■ 

3.3. Attainability. With f|3.6|) we consider, for n = 3 and n > 3 separately, 
the hitting times for l<i^j<n— 1: 

ry := inf {f > : Y(t) G } = inf{i > : Y(t) G } . 

First we look at the case n = 3, i.e., two-dimensional reflected Brownian 
motion and the hitting time t\% of the origin. The directions of reflection ti 
and t2 can be written in terms of angles. Note that the angle £ of the two- 
dimensional wedge 6 is positive and smaller than it, since all the eigenvalues 
of 21 are positive. Let 9\ and 62 with — ir/2 < 9%, 62 < vr/2 be the angles 
between ni and ti, and between H2 and ?2 j respectively, measured so that 9\ 
is positive if and only if ti points towards the corner with local coordinate 
(0,0)'; similarly for 62 . See Figure [Din the end of this paper. 

Paraphrasing the result of Varadhan & Williams for Brownian motion 
reflected on the ta;o-dimensional wedge, we obtain the following result on 
the relationship between the stopping time and the sum 6i + 9j of angles of 
reflection directions, when n — 1 = 2 . 

Lemma 3.1. (Theorem 2.2 of [25]) Suppose that Y(0) = y G 6 \ $° , and 

consider the ratio (3 := {0\ + #2)/£ • 
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The submartingale problem for the reflected Brownian motion on the two- 
dimensional wedge is well-posed for (3 < 2 , whereas it has no solution for 
[3 > 2 . If < j3 < 2 , we have P(ri2 < oo) = 1 ; if, on the other hand, 
(3 < , then we have P(ri2 < oo) = . 

In terms of the reflection vectors ni,ti and ri2,t2 , and with the aid of 
(13. 6j) . we can cast this result as follows; the proof is in subsection 15.41 

Lemma 3.2. Suppose that Y(0) = y G M? \ . If n[q 2 + n 2 qi > , then 
we have P(ti2 < oo) = 1 . If, on the other hand, + n 2 qi < 0, then 

P(ri2 < oo) = . 

We consider the general case n > 3 next. From (|3.6p and Theorem 1.1 
of Williams [26] we obtain the following result, valid for n > 3. 

Lemma 3.3. Suppose that Y(0) = yo G IR™ -1 \ 5° and n > 3, and that the 
so-called skew-symmetry condition 

(3.7) n-qj + = 0; 1 <i < j <n - 1 

holds. Then we have P(r < oo) = , where r := inf{i > : Y(t) G 5° } ■ 

Moreover, the components of the adapted, continuous and non- decreasing 
process L(-) defined in t \3.4\) are identified then as the local times at the 
origin of the one- dimensional component processes: 

2Li(t) =Yi(t)-Yi(0)- [ sgn(y i (s))dy i (s); < t < oo , i = l,...,n. 

Jo 

Remark 3.1. In the planar (two-dimensional) setting of Lemma 13.21 the 
skew-symmetry condition (|3.7p takes a weaker form, that of an inequality. 
In the next section we shall discuss some details of the resulting model, as 
an application of Lemma 13.31 □ 

Lemmata 13.21 and 13.31 lead to the following result, proved in subsection 
14.2.21 on the absence of triple-collisions for a system of n one-dimensional 
Brownian particles interacting through their ranks. Let us introduce a col- 
lection {Qu }i<i k<n °f polyhedral domains in W 1 , such that {Q^ }i<j< n 
is partition W 1 for each fixed k, and {q1 }i<fc<n is partition IR n for each 
fixed i. By analogy with (|3.ip . the interpretation is as follows: 

V = (Vl, ■ ■ ■ , Vn)' G Qj/ means that j/j is ranked kth among j/i, • • • , y n , 

with ties resolved by resorting to the smallest index for the highest rank. 

Proposition 3. For n > 3 , consider the weak solution of the equation 
with diffusion coefficient (ji.ffl) . where a(-) is the diagonal matrix 



(3.8) <t(x) := diag ^ a k^ Q m (x), . . . , ^ a k l (n) (x) ; xG 



k=l 



If the positive constants {ofc;l<£;<n} satisfy the linear growth condition 

Ci n^ ~2 ~2 ~2 ~2 ~2 ~2 

(3.9) cr 2 - o"! = cr 3 - £7 2 = • • • = a n - a n _ x , 
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then holds: there are no triple- collisions among the n particles. 

If n = 3, the weaker condition Tj\ — o\ > a% — <r| is sufficient for the 
absence of triple collisions. 



Remark 3.2. The special structure (I3.8P has been studied in the context of 
Mathematical Finance. Recent work on interacting particle systems by Pal 
& Pitman [19] clarifies the long-range behavior of the spacings between the 
arranged Brownian particles under the equal variance condition: 5"i = • • • = 
a n ; the setting of systems with countably many particle is also studied there, 
and related work from Mathematical Physics on competing tagged particle 
systems is surveyed. The "linear growth" condition (|3.9|) should be seen in 
the light of Figure 5.5, p. 109 in Fernholz [8]. □ 

4. Application 

4.1. Atlas model for an Equity Market. Let us recall the Atlas model 

(n \ n 

Y,9ki Q i,){x{t)) + 1 \dt + J2^ Q < j i)(x(t))dw i (ty, 
k=l / k=l 

for 1 < i < n, < t < oo, (^i(O), • • • ,X n (0))' = x G M n , 

introduced by Fernholz [8] and studied by Banner, Fernholz & Karatzas 
PQ. Here X(-) = (Xi(-), • • • ,X n {-))' represents the vector the logarithms of 
asset capitalizations in an equity market, and we are using the notation of 
Proposition 3. We assume that the constants > and , k = 1, • • • , n 
satisfy the following conditions, which ensure that X(-) is ergodic: 

9i < , g\ + 52 < , • • • , g\ H h g n -i < , g\ H V g n = . 

The dynamics of (|4.ip induce corresponding dynamics for the ranked pro- 
cesses X^(-) > A( 2 )(-) > ••• > A( n )(-) of (|3.ip . These involve the local 
times A k ' e (-) = L x ( k )~ x W(-) for 1 < k < i < n , where L Y '(•) denotes 
the local time at the origin of a continuous semimartingale Y(-) > . An 
increase in A k,i (-) is due to a simultaneous collision of t — k + 1 particles 
in the ranks k through i . In general, when multiple collisions can occur, 
there are (n — l)n/2 such possible local times; all these appear then in the 
dynamics of the ranked processes, as in Banner &: Ghomrasni [2]. 

Let Sk(t) := {i : Xi(t) = X^(t)} be the set of indices of processes 
which are kth ranked, and denote its cardinality by Nk(t) := \Sk{t)\ for 
< t < oo . Banner & Ghomrasni show in Theorem 2.3 of [2] that for any 
n— dimensional continuous semimartingale X(-) = (Xi(-), . . . , X n (-)) , its 
ranked process X^(-) with components X^(t) = X pt ^(t), k = 1, .. .,n is 
(4-2) 

n 1 n fc— 1 

dx (k)(t) = Ehx (k) (t)=x m dX l (t) + 1 ^[ £ dA^)-£dA^(t) 

i=l kK ' j=k+l j=l 
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Here pt := { (pt(l)j • • • ,Pt{n)) } is the random permutation of {1, . . . , n} 
which describes the relation between the indices of X(t) and the ranks of 
such that p t (k) < p t (k + 1) if X^(t) = for < t < oo . 

Let H n be the symmetric group of permutations of { 1, ...,n}. The 
map pt : f2 x [0, oo) i— > II„ is measurable with respect to cr-field generated 
by the adapted continuous process {X(s) ,0 < s < t} , so is predictable. 
Consider the inverse map p^ 1 := (p^~ (1), . . . ,p^ (n)) : 17 x [0, oo) — ► II n , 
also predictable, indicating the rank of Xi(t) in the n— dimensional vector 
X(t) : 

(4.3) X (p _ 1(i)) (t) = X i (t); i = l,...,n, < t < oo . 

Under the assumption of "no triple collisions" (that is, when the only 
non-zero change-of-rank local times are those of the form A fc ' fc+1 (-) , 1 < 
k < n — 1), Fernholz [8] considered the stochastic differential equation of 
the vector of ranked process Xr \ in a general framework; Banner, Fernholz 
& Karatzas [I] obtained a rather complete analysis of the Atlas model (|4.ip . 

In this section we apply the main results of the previous sections to the 
Atlas model. There are some cases of piecewise constant diffusion coeffi- 
cients, which satisfy the conditions in Proposition 1 or 3. Obviously, if {3^} 
are the same, we are in the case of standard Brownian motion. A bit more 
interestingly, if {5^} are linearly growing in the sense of (13. 9p . we can con- 
struct a weak solution to (14. ip with no collision of three or more particles. 

Remark 4.1. On page 2305, the paper pQ contains the erroneous statement 
that the "uniform non-degeneracy of the variance structure and boundedness 
of the drift coefficients" preclude triple collisions. Part of our motivation in 
undertaking the present work was a desire to correct this error. 

4.2. Construction of weak solution. 

4.2.1. Reflected Brownian motion. Let us start by writing the dynamics of 
the sum (total log-capitalization) X(t) := X\ (•) + ■■• + X n {-) as 

n n n 

(4.4) dX{t) = n-fdt + ^2Yj^ kl Q^^ X ^ dWi ^ =njdt + ^2^kdB k {t), 

i=l k=l " k=l 

where B(-) := {(B 1 (t),--- ,B n (t))',0 < t < oo} is given by B k {t) := 
ELi/n Uj^fs))^^) for 1 < k < n, < t < oo . By the F. Knight 

theorem (e.g., Chapter 3 in Karatzas & Shreve [E]), this process B(-) is 
n— dimensional Brownian motion started at the origin. 

Next, let h and £ be the (n-1) X 1 vector and the (n — 1) x n triangular 
matrix with entries 

hi V v ?2 " ?3 

h •■= {91-92,--- ,9n-l-9n) , ^ : = 

\ <?n-l -a n J 
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where the elements in the lower-triangular part and the upper-triangular 
part, except the first diagonal above the main diagonal, are zeros. Then 
the process {ht + ££?(£), < t < 00} is an (n — 1)— dimensional Brownian 
motion starting at the origin of 
matrix 



(4.5) 



21 



EE' 



V 



-(To 



n 1 , with constant drift h and covariance 

\ 



al + al 



'n-l 



'n-l 



J 



Now we construct as in Section I3T21 an (n— 1)— dimensional process Z{-) :- 

n-l 



{(Zi(t), ■■■ , Z n -i(t))' , < t < 00 } on R™- 1 by 

Zk(t) ■= {gk - 9k+i)t + &kBk(t) - a k+ iB k+ i(t) 

+ A k ' k+1 {t) - ^ (A* _1 ' fc (t) + A fc+1 ' fc+2 (i)) ; < t < 



(4.6) 



:3C 



for k = 1, • • • — Here A k,k+1 (-) is a continuous, adapted and non- 
decreasing process with A fe ' fc+1 (0) = and / °° Z k (t) dA k ' k+1 (t) = almost 
surely. Setting A 0,1 (-) = A n,n+1 (-) = 0, we write in matrix form: 



Here A(- 



(4.7) 



Z(t) 



I-Q 



/it + E J B(t)+9 c tA(t); 0<t<oo, 

, A k ~ 1,k (-))' and the reflection matrix IK 

/ 1 -V2 \ 
-1/2 1 



••• -1/2 
-1/2 1 



J-jQ is 



If the process X(-) has no "triple collisions", then from (|4.2p we get 
dX (k) (t) = ^! {XiW=X(fc)W} dA,(t)+-(dA^ +1 (t)- ( iA fc - 1 ' fe (t)) , < t < 



CXD 



i=l 



Substituting (|4.ip into this equation and subtracting, we obtain that 

(4.8) X (k) (t) - X {k+1) (t) = Z k (t); l<k<n-l, 0<t<oo, 

and that A fc ' fc+1 (-) is the local time at the origin of the one-dimensional 
process Z k {-) > for k = 1, . . . n— 1 . In general, the process A(-) may have 
triple (or higher- multiplicity) collisions, so that we have additional terms in 



(4.9) X (k) (t) - X (k+l) {t) = Z k {t) + Ck(t) , l<fc<n-l,0<i<oo. 
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The contribution £(•) := (Ci( - )> • • • > Cn-i( - )) from triple or higher-multiplicity 
collisions can be written for 1 < k < n — 1 , < i < oo as £fc(0) = and 

n n k—2 

£=3 jf=fc+2 j=l 

n n k—1 

£=3 j=fc+3 i=l 

Remark 4.2. Note that £(•) consists of (random) linear combinations of local 
times from collisions of three or more particles. It is flat, unless there are 
triple collisions: i.e., f Q = 0, where <5 := {s > : X,i{t) = 

Xj(t) = X k (t) for some l<i<j'<A;<n}. We use this fact with Lemma 
14.11 in the next subsection. □ 

4.2.2. Proof of Proposition 3. Under the assumption of Proposition 3, we 
can apply Lemma 13.31 to obtain 

(4.10) ¥(Zi(t) = Zj{t) = , 3t > , , 1 < j < n) = ; 

see subsection 15.51 Thus Z(-) is a special case of multi-dimensional re- 
flected Brownian motion for which each continuous, non-decreasing process 
A k,k+1 (-) is exactly the local time at the origin of Zfc(-) . 

Now let us state the following lemma to examine the local times from 
collisions of three or more particles. Its proof is in subsection 15.61 

Lemma 4.1. Let a(-) = {a(t) ; <t < oo } be a non-negative continuous 
function with decomposition a(t) = f3(t)+j(t) , where /?(•) is strictly positive 
and continuous, and 7(-) is of finite-variation and flat off {t > : a(t) = 
0} , i.e., J Q l/ a m > o}d7(f) = 0. Assume 7(0) = and a(0) = (3(0) > 0; 
then, 7(i) = and a(t) = (3(t) for all < t < 00 . 

Under the assumption of Proposition 3, applying the above Lemma 14.11 
with (j49]), flUUD and a(-) = X (k) (;u) - X (k+2) (; u>) , /?(•) = Z k {-,u) + 
Z k+ i(-,u) and 7Q = Cfe(-,^) + Cfc+l(-)^) for u G CI , we obtain a(-) =/?(•) : 

(4.11) X {k) (-) - X {k+2) (-) = Z k (-) + Z k+1 ( : ) , k = l,...,n-2. 
Combining (|4.11j) with (14.10j) . we obtain X^(-) - X( fc+2 )(-) > or 

F(X (k) (t) = X (k+1) (t) = X (k+2) (t) , 3t > , 3k , 1 < k < n - 2) = . 

Therefore, there are "no triple collisions" under the assumption of Proposi- 
tion 3, whose proof is now complete. □ 

4.2.3. Recovery. In conclusion, we recover the n-dimensional ranked process 
Xr\ of X by considering a linear transformation. Specifically, we construct 
the n— dimensional "ranked" process 

tf (0 (t) := (*(!)(*),••• ,*(*)(*)); 0<t<oo 
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from the sum X(t), < t < oo denned in (j4.4H and the reflected Brownian 
motion Z(-), so that the differences (gaps) satisfy 

(4.12) *( fc )(t) -*(*+!) (*) = fc = l,...n-l, 
and the sum satisfies 

n 

(4.13) £tf (fc) (i) = £(*); 0<t<oo. 

k=l 

In particular, each component of ^f^(t) is uniquely determined by 
/ * (1) (t) \ / + Z n _i(t) + (n - 2)Z n _ 2 (t) + • •• + (n - l)Zi(t) \ 



*(2)(*) 

v *(«)(*) y 



£(t) + Z n _i(-) + (n- 2)Z n _ 2 (t) + Zx(t) 



V 3E(t) - (n- l)Z n _i(t) - (n- 2)Z„_ 2 (*) Z x (t) / 

for < i < oo. Under the assumption of Proposition 3, we obtain (|4.10p 
and hence with (|4.12p we arrive, in the same way as discussed in f|3.3[) . at 

P(* (fc) (i) = ¥( fc+ i)(t) = *(fc +2 )(t), 3t > 0, 1 < 3& < n - 2) = 0. 

Thus, the ranked process {Xr.\(t),0 < t < oo } of the original process X(-) 
without collision of three or more particles, and the ranked process ^(.)(-) 
defined in the above, are equivalent, since both of them have the same sum 
(|4,13p and the same non- negative difference processes Z(-) identified in (|4.8p 
and (|4,12p . We may thus view as the weak solution to the SDE for the 

ranked process Xr.\(-) . Finally, we define := (*].(•)> • • • ;^n( - )) where 
^i(-) = for i = 1, . . . , n , and (i) is defined in (|4.3p . Then, 

^(•) is the weak solution of SDE (|4.ip . This construction of solution leads 
us to the invariance properties of the Atlas model given in pQ and |13j . 

5. Appendix 

5.1. Proof of Lemma 12.11 From the assumption xo € M n \ Z , where the 
zero set Z is defined in (|2.10p . it follows that s(0) = s(X(A )) > and 
there exists an integer niQ such that 1 < s(0) < mo ■ Recall that with 
R(X(A.)) = D(-) and s(X(A.)) = s(^) we obtained (I2TT41 . namely: 

s (t) = s(0) + [ *( u \~ du + B(t) ; < t < oo. 
Jo 2s (u) 

• Let us consider first the case d := essup sup 0<4<00 D(t) < 2 for (|2.16p . 
Define two continuous functions b\(x) := (0 — 1) / (2x) and ^(^O :=f> / (4x) 
for x G (0, oo) . If d < 2 , then &i(-) < &2( - ) in (0, oo) . For each integer m > 
mo , there exists a non- increasing Lipschitz continuous function / m (-) := 
(foi (•) + 62 (•)) / ^ with Lipschitz coefficient K m := max^g^-i^l&^x)! , such 
that 6i(.) < fm{-) < 6a(-) m [m-\m] . 

Define an auxiliary Bessel^ process r(-) of dimension (0 + 2) / 2 (< 2) : 

x(t):=s(0)+ b 2 (x(u))du + B(t); < t < 00 . 
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Consider also the increasing sequence of stopping times 

(5.1) T m : = inf{i > : max[s(i),t(f)] > m or min[s(i), t(i)] < mT 1 } 

for rriQ < m < oo , and r := inf{i > : t(t) = } . From the property of the 
Bessel process with dimension strictly less than 2, the process t(-) attains 

the origin within finite time; : — lim m >00 T m < t < oo holds a.s. 

Now take a strictly decreasing sequence {a n }^ =0 C (0, 1] with ao = 1 , 
linin^oo a n = and J^ a a ^u~ 2 du = n for every n > 1. For each 

n > 1 , there exists a continuous function p n (-) on R with support in 
(a n ,a n _i), so that < p n (x) < 2(nx 2 )~ l holds for every x > and 
/(a n ,a„_i) z 9 ™^) dx = 1 - Then the function ip n (x) := ( f£ p n (u) d u) d y; 
i £ I is even and twice continuous differentiable with |^(x)| < 1 and 
lim^oo ip n (x) = \x\ for Define tp n (') := ip n (-) l(o )00 )(-) • 

By combining the properties of tp n (-) , &i (•) , &2 (•) and f m { ) , we see that 
the difference A(-) := s(-) — r(-) is a continuous process with 

M^(t)) < ! y/ n (A(«))(6i(B(«)) - 6a(t(u))) 

■/ o 

<[ ip' n {&{u)){f m (s{u))-f m {t{u)))du<K m ! ^(A(n))(s(n)-r( U )) + d U 

JO ■/ 

<K m [ (A(u)) + du; 0<t<r m . 
J o 

Letting n — > oo we obtain (A(i)) + < X m JjJ (A(u)) + dit for < i < r m . 
From the Gronwall inequality and he sample-path continuity of s(-) , t(-) in 
[0, oo) , we obtain A(-) = s(-) — r(-) < on [0, r m ] for m > tuq and 

(5.2) s(t*) = lim s(i) < lim x(t) = r(r*) and max[s(r*),r(r*)] < oo , 

t— >T* t— »T* 

almost surely. On the other hand, from the definition of {r m } we obtain 

= t(r*) > s(t*) , thus s(r*) = and s(t) < x(t) for < t < r* , a.s., so 
for i) = essup sup 0<i<oo i)(t) < 2 we conclude 

Q X0 ( s(X(t)) = for some t > 0) = ( a(t) = for some t > ) = 1 . 

By the strong Markov property of the process X(-) under Q , we obtain 

1 = Qx (s(X(t)) = , inf. many t>0)= Q xo ( s(t) = , inf. many t > 0) . 

This gives ()2. 16j) of Lemma [2.1i Moreover, by the formula of the first hitting- 
time probability density function for the Bessel process with dimension d 
in Elworthy et al. [7j and Going-Jaeschke & Yor [10] we obtain 

Q X0 (s(t) = , for some t G (0, T] ) 

> Q xo {t(t)=0, for some t£ (0, T] ) = 1 - k(T; s(x ),d) , 

where the tail probability distribution function «(•;•,•) is defined in (I2.18|) . 
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• We consider next the case of c) := essinf info<t<oo 0(i) > 2. Define 

bz(x) := (0 — 1) / (2x) and b^{x) := Z> / (4a;) for x G (0, oo) . Following 

a course similar to the previous case, using 6 3 (-) , b±{-) and defining a non- 
increasing Lipschitz continuous function g m { ) := (6 3 (-) + &4(-)) / 2 with Lip- 
schitz coefficient L m := maxjgu-i^i | &3 (ar) ] (rather than using &]_(•) , &20 > 
/ m (-) and if m ), we obtain the reverse inequality q(-) < s(-) in [0, r m ] a.s. 
Here q(-) is the Bessel process in dimension (9 + 2) / 2 (> 2), namely 

q(i) =s(0) + / 6 4 (q(u))citi + ; < t < oo , 

and the stopping times {r m } are defined as in (|5.ip but with t(-) replaced 
by q(-) . By a well-known property for Bessel processes of dimension at least 
2, the process q(-) never attains the origin: i.e., q(-) > on [0, oo) , a.s. 

If r* := linim^oo r m < oo , then by analogy with (|5.2p we obtain s(%) > 
q(T*) > and max[s(r*) , q(T*)] < oo a.s., and from the construction of 
{r m } a contradiction follows: = s(r*) > 0. Therefore, Q Xo (s(i) > 
for < t < oo ) = 1 . This gives (|2"7T5|) of Lemma [2J] for D > 2 . □ 



5.2. Proof of Propositions 1, 2. Proposition 1 and the first half of Propo- 
sition 2 are direct consequences of Lemma \2. II and of the reasoning developed 
in subsection 12.21 Note that (&}(t) > c^t ,t > in this uniformly nonde- 
generate case. We obtain (|2.22p . because Q XQ (s(X(t)) = for some t E 



[o,r]) > q X0 (s(u 

ability measure IF 
semimartingale with decomposition 



= 0, for some u S [0, cqT]) . Under the original prob- 
Xo , because of the drift /u(-) , the process s(X(-)) is a 



ds(X(t)) 



(R(x) - l)Q(s) x'DD'njx) 



\ 2s(x) s(x) I 

h(X(t))dt + dO(t); < t < oo 

§(■ 



dt + dQ(t) 



x=X(t) 



where /i(-) , ©(•) are obtained from h(-) , &(■) in ()2.12p . upon replacing 
R(-) in dZTB by J2(.) in (I2T23P and W(-) in {£2} by W(-) . The com- 
parison with Bessel processes is then repeated in a similar manner. When 
sup^Rn^ R(x) < 2 , we get (^2ll) and (|2T25|) . 



5.3. Example in Remark 12.31 With some computations we obtain the 
following simplification of the effective dimension given in (I2.29|) : 



EB A (x) 



2 + ^ 



\x\\ 2 - 4ai 2 (x) • xix 2 l^ 1+ u^i_ 
-4a 23 (x) • x 2 x 3 1tz 2+u tz 2 _ 
-4a 3 i(x) • x 3 xil^3 + u^ 3 _ 

x'A{x)x 



for x e R 3 \ {0} 



and 
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4a i2 (x) • [(xi - x 2 ) 2 + 2(x 2 - x 3 )(x 3 - xi)]LK 1+ ufti_ 
+ 4a 23 (x) • \{x 2 - x 3 ) 2 + 2(x 3 - xi)(xi - x 2 )]l^ 2+ u7J 2 - 
_ 9 \ + 4a 3 i(x) • [(x 3 - xi) 2 + 2(x 2 - x 3 )(xi - x 2 )}liz 3+ un 3 - 
{X) ~ + x'DD'A{x)DD'x 

for x G R 3 \Z , where #(•) is defined in (EHfflj) and Z(-) is denned in (l2~TU1) . 
Under the specification (|2.28j) . we verify ED(-) > 2 and R(-) > 2, since the 
denominators of the fractions on the right-hand sides are positive quadratic 
forms and their numerators can be written as 

||x|| 2 — 4ai2(x)xiX2 = (1 — 4a 2 2 )x 2 + xf + (xi — 2a\ 2 x 2 ) 2 > ; x G 7Zi + U Hi- 
4ai 2 (x)[(xi - x 2 ) 2 + 2(x 2 - x 3 )(x 3 - x x )\ = 4ai 2 (x)fi(x) > ; x G TZ 1+ U Ki- 
with similar formulas for x G 7£j+ U 7£j_ , i = 2, 3 . 

5.4. Proof of Lemma 3.2. We recall the special geometric structure of 
orthogonality = and ||nj|| = 1, and observe that 

(5.3) (vt'Q + Q'Vt). ^ ^ nfo + n$q, ^ ; V 

Note that if n = 3 , i.e. n — 1 = 2, then n^qj = ||cjj|| sgn(— #.,■) sin(£) for 
1 < « / j < 2, where sgn(x) := 1{ X>0 } - 1{ X< 0} ■ The length ||q 2 || of q 2 
determines the angle 9 2 and vice versa, i.e., 

~ IRill ~ \Qj\- 

With this observation and < £ < 7r, sin(£) > 0, we obtain 

n-q j +n^q i = sin(0(||qj||sgn(-^) + ||qi||sgn(-6» i )) ~ 

<=> p = {e i + e j )/( ^o; i<i^j<2. 



Thus, we apply Lemma 13.11 and obtain Lemma 

5.5. Coefficient Structure, and Proof of (14.101) . Next, we consider the 
case of linearly growing variance coefficients defined in (|3.9p , and recall the 
tri-diagonal matrices 21 = ££' as in (|4.5p and 9\ as in (|4.7p . Consider the 
(n — l)-dimensional reflected Brownian motion Y(-) defined in (|3.4p with 
S = £ and this above 9^. Such a pair (S,9\) satisfies 



(5.4) 2£> - 09 - 2)0 - 221 =0; l<i,i<n-l, 

V / 

where 2) is the diagonal matrix with the same diagonal elements as 21 of 
(13. 5p . and is the ((n — 1) x (n — 1)) matrix whose first-diagonal elements 
above and below the main diagonal are all 1/2 and other elements are zeros 
as in (14. 5j) . In fact, it suffices to consider j = i + 1, i = 2, ... ,n — 1, for 
which the equalities (|5.4h are 

- (5? + 5 2 +1 ) - (2f 2 _! + 5?) + 45 2 = , 
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or equivalently (j3.9f) : Sf - a 2 _ x = a 2 +1 - <r 2 for 2 < i < n — 1 . 

Moreover, the equalities (|5.4f) are equivalent to + Q'y\)ij = in 

(pOj) . In fact, from (J33]) with S 1 / 2 = GT 1 we compute 

= S)- 1 /2(/ _ £2)2)1/2 _ ^-1/2^-1/2 
9t'£2 + Q'm = 21- D _1 / 2 £2D 1 / 2 - D^jQST 1 / 2 - 2T)- 1/2 ^D- 1/2 , 

and multiply both from the left and the right by the diagonal matrix D 1 / 2 
whose diagonal elements are all positive: 

(5.5) D 1/2 (^'0 + Q!%& /2 = 22) - O.D - DQ - 221. 

The equality in the relation (|5.4p is equivalent to the so-called skew- symmetry 
condition 9VH + £}'9T = introduced and studied by Harrison & Williams 
in [12], [26]. It follows from (pT3j) . ([53D and (j53|) that the reflected Brow- 
nian motion Z(-) defined in (|4.6f) . under the assumption of Proposition 3, 
is such that any two dimensional process (Zi,Zj) never attains the corner 
(0, 0)' for 1 < i < j < n — 1 , i.e., \4- Mty holds. Using this fact, we construct 
a weak solution to (|4,ip from the reflected Brownian motion. This final step 
is explained as an application in the last part of subsection 14.2.21 

5.6. Proof of Lemma 14.11 We fix an arbitrary T £ [0,oo) . Since /?(•) is 
strictly positive, we cannot have simultaneously a(t) = (5{t) +j(t) = and 
7(t) > 0. The continuous function (3{-) attains its minimum on [0, T] , so 

{t e [0,T] : a(t) = 0} = {t G [0,T] : a(t) =0, j(t) < 0} 

( 5 - 6 ^ C {t e [0,T\ : 7(t) < - min /3(s) < 0}. 

0<s<T 

Let us define to := inf{i G [0,T] : a(t) = 0},with to = 00 if the set is 
empty. If to = 00 > then a(t) > for < t < oo ; thus, it follows from 
the assumptions 7(0) = and f Q l{a(t)>o}d^/(t) = for < T < 00 that 
7(-) = . On the other hand, if to < 00 , then it follows from the same 
argument as in (|5.6f) that 7(to) < — mino< s <t (3(s) < . This is impossible, 
however, since a(s) > for < s < to by the definition of to and hence the 
continuous function 7(-) is flat on [0, to), i.e., = 7(0) = 7(to~ ) = 7(^0) • 
Thus, to = 00 and 7(-) = . Therefore, the conclusions of Lemma I4TT1 hold. 
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Figure 1 . Directions of reflection: 6\ + 9 2 < 
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Figure 2. Polyhedral regions in Remark 12.31 
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